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In this work we report a new route to chaos from a resonance torus in a piecewise smooth non- 
invertible map of the plane into itself. The closed invariant curve defining the resonance torus is 
formed by the union of unstable manifolds of saddle cycle and the points of stable cycle and saddle 
cycle. We have found that a cusp torus cannot develop before the onset of chaos, though the loop 
torus appears. The destruction of the two-dimensional torus occurs through homoclinic bifurcation 
in the presence of an infinite number of loops on the invariant curve. We show that owing to the non- 
invertible nature of the map, the structure of the basin of attraction changes from simply connected 
to a nonconnected one. We also describe how the mechanism of transition to chaos differs from the 
scenario of appearance of chaos in invertible maps as well as in smooth non-invertible maps. 

PACS numbers: 05.45.Gg, 05.45.Pq 



I. INTRODUCTION 

Piecewise smooth (PWS) maps have received much re- 
search attention in the recent times because of their ap- 
plicability in a large number of practical systems [3, [2, 0] • 
This includes switching circuits [J], impact oscillators 
[sl EJ , a variety of micro and macro economic systems @ , 
cardiac dynamics [8[ and many other systems. For such 
maps a particular kind of bifurcation may occur, which 
is different from what can occur in smooth maps. As a 
control parameter changes, the fixed point may move in 
the phase space and may collide with the borderline be- 
tween two smooth regions, resulting in an abrupt jump 
of the Floquet multipliers from the inside of the unit cir- 
cle to the outside of it in the complex plane, leading to 
a special class of nonlinear phenomena known as border 
collision bifurcation (BCB) 0, i, [TO, EL El . 

The BCB consists of many atypical phenomena such as 
direct transition from a period- 1 attractor to a chaotic at- 
tractor [1] , m ultiple attractor bifurcation [l3[ and danger- 
ous BCB [14]. Recently in a series of papers [l5l [l6[Tl7t 
a new type of BCB has been reported, where a stable 
fixed point changes into an unstable focus along with 
an attracting closed invariant curve. This closed invari- 
ant curve forms a torus associated with quasiperiodic 
or phase- locked periodic dynamics [Hj]. Moreover, the 
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torus breakdown route to chaos in PWS invertible map 
has been characterized in 16] which is distinct from that 
of Afraimovich-Shilnikov [18|. These mechanisms of the 
creation and destruction of tori were studied in the con- 
text of two-dimensional smooth invertible maps [l9[ and 
PWS invertible maps T$\. 

In another line of work [2(J [U, [22[ , the occurrence of 
chaos via torus destruction has been studied numerically 
and experimentally for non-invertible smooth maps. Sev- 
eral routes to chaos were reported. The basic theory was 
proposed by Maistrenko et al. [20[ and Lorenz [21|. In 
transition from a smooth invariant curve (IC) to a chaotic 
attractor, the process begins with the IC developing seg- 
ments with increasingly high curvature. At a certain 
parameter value cusps develop on the IC, which subse- 
quently change into loops. With further change of the 
parameter, the IC can be destroyed in accordance with 
scenarios analogous to those of the invertible case. After 
that the chaotic attractor comes into existence. Cusps 
appear on IC when the degree of a map is greater than 
one and the IC intersects the critical curve so that the 
tangent of the IC at the point of intersection coincides 
with the eigenvector corresponding to zero eigenvalue or 
the critical curve contains cusp point (a point on a curve 
at which tangents of each branch of the curve coincide). 
Critical curve is the two dimensional extension of the 
local extrema (critical point) of a one dimensional differ- 
entiable map. It is the locus of all those points at which 
Jacobian determinant vanishes. For a non-differentiable 
piecewise linear map like the tent map, the point of non- 
differentiability plays the role of critical point. But for 
a two dimensional piecewise smooth map, according to 
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Mira et al. [23] the critical curve LC is the image of 
the borderline LC-i of the map and is characterized by 
the set of all points which have at least two coincident 
rank-1 preimages. This curve may contain more than one 
segment. It separates the phase plane into open regions 
where all points of a region have the same number of 
rank-1 preimages. 

The two dimensional PWS normal form maps 0, EH 
investigated so far are invertible. In the present paper we 
are interested in the dynamical behavior of non-invertible 
PWS maps. Considering a general piecewise linear 2D 
map expressed in the normal form, we derive the pa- 
rameter regions where the map is non-invertible. We il- 
lustrate how chaos appears in PWS non-invertible maps 
and demonstrate the difference from the corresponding 
mechanisms of transition for smooth invertible [19] and 
non-invertible maps [gO, Hl| . We investigate the bifurca- 
tion that takes place as one moves from the inside of a 
resonance tongue to the outside in the parameter space. 
Inside a resonance tongue, a closed invariant curve is 
formed by the saddle-node connections of a pair of cycles 
(a saddle and a node). This curve forms the mode- locked 
torus. We show that the cusp torus does not develop 
but the loop torus [20[ is created before the destruction 
of mode-locked torus as the parameter is varied contin- 
uously. The mode-locked torus is destroyed through a 
homoclinic bifurcation. We observe that the interaction 
between stable manifold and the critical curve leads to 
a qualitative change in the basin of attraction, from a 
simply connected to a nonconnected one. We also show 
that the unstable manifolds can have structurally sta- 
ble self-intersections while the stable manifolds can in- 
tersect itself when the map becomes structurally unsta- 
ble. There is a coexistence of stable node and a chaotic 
attr actor. We show that the size of basin of attraction of 
the stable cycle (node) diminishes to zero and after this 
the stable node undergoes a border collision fold bifur- 
cation. We also illustrate the important role played by 
non-invertibility of the map in the mechanism of transi- 
tion to chaos. 



II. CONSIDERED SYSTEM 

For a piecewise smooth map whose leading-order Tay- 
lor term close to the border is linear, the dynamics in the 
neighborhood of the border can be expressed by a normal 
form map 0, H defined by 



F : 



where 
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FIG. 1: Bifurcation diagram obtained by varying the parame- 
ter fi from —0.01 to 0.01. This diagram illustrates the birth of 
a quasiperiodic orbit after the border collision, directly from 
a stable fixed point for 5l — —0.3, 5r = 1.4, tl = —0.1 and 
tr = -0.685. 



The phase plane R 2 is divided into two compartments 
L := {(x,y) e M 2 : x < 0} and R := {(x,y) G M 2 : x > 
0}. tl and Sl are the trace and determinant respectively 
of the Jacobian matrix Jl in the left half-plane L, and 
tr and Sr are the trace and determinant respectively of 
the Jacobian matrix Jr in the right half-plane i?, and \i 
is a parameter. Thus 



Jl 



TL 1 

-S L 



and Jr 



TR 1 

Sr 



Most of the theories of BCB discussed so far are on 
dissipative dynamical systems (i.e., \5l\ < 1, \Sr\ < 1). 
In contrast, we consider a situation where the map is 
contracting in the left side and expanding in the right side 
of the border, which happens when \Sl\ < 1, \Sr\ > 1. 
We choose tx, tr satisfying the conditions 

-(1+5 L ) <t l < (1+S L ), -2v^ < tr < 2V^- (2) 

The above conditions ensure that the fixed point for 
\i < is attracting and for \i > it is a spiral repellor. 
According to the notations in [ll| , we can write the stable 
fixed point by A and unstable fixed point by a in L, and 
similarly by B and b respectively in R. The fixed point 
in L is 



A or a = ( : 



' 1 - T L + Sl ' 1 - T L + Sl 
and the fixed point in R is 



B or b = ( 



1-tr + 5r' 1-t r + 5r 



) 



A or a exists if /i/(l — tl + ^l) < 0, otherwise a virtual 
fixed point is located in R and is denoted by A or a. 
Similarly B or b exists if /i/(l — tr + Sr) > 0, otherwise 
a virtual fixed point is located in the L side and is denoted 
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FIG. 2: Regions of the periodic tongues in the (tl,tr) pa- 
rameter plane for Sl = —0.3, Sr = 1.4 and \i — 0.05. Distinct 
colors of the tongues denote different oscillatory motions. 

by B or b. The stability of the fixed points is given by 
the eigenvalues 

Ai,A 2 = i(T±Vr 2 -4£). 

As the parameter \i hits the border at \i = 0, a stable 
fixed point A on the L side becomes an unstable focus 
b on the R side. When the outward spiraling orbit in 
the R side touches the boundary of the L side then it 
is attracted by the virtual fixed point A located in R. 
So the outward motion is arrested. This gives rise to the 
existence of an invariant closed curve on which stationary 
state must lie. Fig. [T] displays the bifurcation diagram 
varying the parameter fi from a negative to a positive 
value. This diagram depicts an abrupt transition from 
periodic to quasiperiodic behavior at a border-collision 
bifurcation. 

We plot the different resonance tongues in (tl,tr) 
parameter plane for a positive value of fi (see Fig. [2]). 
The main resonance tongues are denoted with the cor- 
responding rotation numbers. Inspection of this figure 
reveals the presence of different periodic tongues which 
are densely spread in the parameter space. Yang and 
Hao [24[ described the lens-like resonance tongues for 
a one-dimensional piecewise linear system, which were 
later generalized for two and three dimensional piecewise 
smooth systems in 0, [25|, [26| . The tongues are bounded 
by the curves where the stable and unstable cycles merge 
and disappear through BCB. A tongue may contain more 
than one lens. 



III. PROPERTIES OF THE CONSIDERED 
NON-INVERTIBLE MAP 

Proposition 1. The two dimensional piecewise smooth 
normal form map (Qp is non-invertible, if either Sl and 
Sr are of opposite sign, or any one of Sl and Sr is zero. 



Proof. Since our map is continuous in R 2 , so we have to 
derive the conditions from one to one property. Let us 
suppose 

X 1= X 2 = and*!, X 2 eR 2 . 

There are three possibilities. 
Case I: X ly X 2 G L. 

NowFi(li) =F 1 (X 2 ) 
=> r L xi + yi + fj, = r L x 2 + 2/2 + M an d - S L xi = -S L x 2 . 

xi = x 2 (if S L ^ 0) and y 1 = y 2 . 

50 F\ is one-one if Sl 7^ 0. 
Case II: X u X 2 G R. 

By a similar argument as in Case I, we conclude that 
F 2 is one-one if Sr ^ 0. 

Case III: X x G L, X 2 G R. 

Fi(Xi)=F 2 (X 2 ) 

=> t l xi + 2/1 + /i = t r x 2 + y 2 + /i and - S L xi = -S R x 2 . 
^x 1 = (S R /S L )x 2 (if S L ^0). 

(a) If Sl and Sr are of the same sign and are non-zero. 
Since x\ < and x 2 > 0, the relation x\ — (Sr/Sl)x 2 

cannot hold. Therefore there can exist no points X\ G L, 
X 2 G R such that Fi(Xi) = F 2 (X 2 ). Hence the map is 
invert ible. 

(b) If Sl and Sr are of opposite sign. 

The condition x\ = (Sr/Sl)x 2 may hold. 

Then Fi(Xi) = F 2 (X 2 ) => X x ^ X 2 because x x < 
and x 2 > 0. Hence the map is non-invert ible. Combining 
the above three cases the map is non-invertible if one of 
the following relations is satisfied 

(i) Sl and Sr are of opposite sign; 

(ii) any one of Sl and Sr is zero. □ 

To consider a typical non-invertible case, we assume 

51 = —0.3 and Sr = 1.4. Each point in R 2 is mapped 
under (pQ) to a point in R 2 with ^/-coordinate less than 
or equal to zero. There are two inverses for (x, y) in R 2 
with y < 0. These are defined as 

F-\x,y)=( Jr 1 ^ ) and 
1 v ' yj \x + (T L /S L )y- v J 

2 v ' yj \x + (T R /S R )y- fi J ' 

where (x, y) G R 2 and y < 0. 

Hence a point in R 2 has two distinct rank-1 preimages 
if y < 0, no preimage if y > and two coincident preim- 
ages if y = 0. Let Zo be the region with no preimage (the 
upper half plane) and Z 2 be the region with two different 
preimages (the lower half plane). Thus for our choice of 
parameters, the map is of Zo — Z 2 type [23[. The locus 
of the points having two merging rank-1 preimages is the 
x-axis, called LC. Here LC is the critical curve. The 
preimages of each point in LC lie on the 2/-axis, named 
as LC-i. So the non-invertible map folds the phase plane 
along LC. The fixed points, periodic orbits and attractor 
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must lie inside Z 2 , because every fixed point and every 
point of a periodic orbit have an inverse (namely itself). 
A point on an attractor also has a preimage lying on that 
attractor. 



IV. STRUCTURE OF THE STABLE AND 
UNSTABLE MANIFOLDS 

Let / : R 2 — > R 2 be a non-invertible smooth map 
and x be a hyperbolic period-^ point of /. We de- 
fine the stable manifold of x as W s (x) = {y e R 2 : 
f nq (y) — » x as n —> 00 } an d the unstable manifold of 
x as W u (x) = {y e R 2 : f~ nq (y) — > x as n — > 00}. 
These are defined in the neighborhood of x, and hence 
so are actually the local stable and unstable manifolds 
of x [27]. The global unstable manifold is formed by 
the union of forward images of the local unstable mani- 
fold. It is connected, since the local unstable manifold is 
connected and the image of a connected set under con- 
tinuous map is connected. The global unstable manifold 
is forward invariant but may not be backward invariant 
(some of its points have multiple preimages) and so self- 
intersections are allowed. The global stable manifold is 
formed by the union of all preimages of any rank of the 
local stable manifold. Since the map / is non-invertible, 
the global stable manifold is backward invariant, but it 
may be strictly mapped into itself (some of its points 
may have no preimage). So it may be disconnected. For 
non-invertible maps the use of the term "manifold" is an 
abuse of the terminology; the local stable and unstable 
manifolds are guaranteed to be true manifolds, but the 
global manifolds are not (27], [28[. Still for our map we 
denote the global stable and unstable manifolds as the 
stable and unstable manifolds respectively. 

Since our map F is of Z$ — Z2 type and the local un- 
stable manifold originates at the saddle point x, it must 
lie in the Z2 region. The unstable manifold is grown 
by iterating the local unstable manifold forward in time, 
every point computed will have a preimage, namely the 
previous point iterated to arrive at that point and hence 
the unstable manifold will never leave the Z2 region (see 
Fig. 02(a)). The stable manifold is the union of all preim- 
ages of local stable manifold, connected or disjoint, so it 
can lie inside the Zq region also (see Fig. 02(a)). 



V. TRANSITION TO CHAOS IN THE 
PRESENCE OF HOMOCLINIC STRUCTURE 

Let us now analyze the bifurcational transition from 
mode-locking to chaos. To illustrate this scenario, we 
vary the bifurcation parameter tr from the inside of the 
1 : 4 tongue to the outside of it along the direction D 
(see Fig. [2j). The reason for this choice is that, the tran- 
sition to chaos maintains the loop like structure as de- 
scribed by Maistrenko et al. [2(j which we want to ex- 
plain. Fig. [H(a) displays the bifurcation diagram and 




-0.04 x -0.01(b) 



FIG. 3: (a) Phase portrait of the attractor at tr = 0.25 
with t l -0.3, S l -0.3, Sr 1.4 and fi 0.05. (b) 
Enlargement of a portion of the IC marked by the square in 
(a). Points of stable cycles are marked with solid black circles, 
points of saddle cycles with open circles and the spiral repellor 
is marked with solid green circle. 



Fig. 21(b) depicts the Lyapunov exponents corresponding 
to the bifurcation diagram. The largest Lyapunov ex- 
ponent becomes positive at tr « 0.5, which implies the 
existence of a chaotic attractor. The second Lyapunov 
exponent remains negative over the considered param- 
eter range. Initially the map (pQ) has a period-4 stable 
cycle and a period-4 saddle cycle. As the parameter tr 
increases, between ~ 0.5 and ~ 0.57414 the stable 
period-4 cycle coexists with a chaotic attractor. So we 
observe hard hysteretic transitions from periodic mode to 
chaotic mode and vice versa, at the points and re- 
spectively. At one of the points of the period-4 stable 
cycle and a point of the period-4 saddle cycle hit the bor- 
der together. Then the period-4 stable and saddle cycles 
disappear through border collision fold bifurcation. 

Let us now characterize the mechanism of occurrence 
of chaos in detail. The phase portrait at tr = 0.25 is 
shown in Fig. [31(a). The period-4 attracting cycle has 
real eigenvalues. The period-4 saddle cycle always has 
one positive eigenvalue Ai > 1 and another eigenvalue 
— 1 < A2 < 0. Since the stable eigenvalue of the sad- 
dle cycle is negative, so our map is orientation revers- 
ing and hence upon successive iterations, points on the 
stable manifold oscillate about the saddle point as they 
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FIG. 4: (a) Bifurcation diagram illustrating the birth of a 
chaotic orbit via border- collision saddle-focus bifurcation and 
(b) the corresponding Lyapunov exponents. The bifurcation 
parameter tr changes from 0.25 to 0.65, for tl — —0.3, 5l — 
—0.3, Sr = 1.4 and [i — 0.05. The period-4 orbit coexists 
with the chaotic orbit in the region t r < tr < r R , where 
t£ « 0.5 and r R « 0.57414. 



approaches to it. The closed invariant curve is formed by 
the union of the unstable manifolds of period-4 flip saddle 
cycle and the points of the period-4 regular attractor and 
flip saddle cycles. It is non-smooth, since the unstable 
manifolds are as smooth as the original map and our map 
is piecewise smooth. This curve forms the mode-locked 
torus. Fig. [3](a) shows that the unstable manifolds are 
self-intersecting. As our map is an endomorphism, the 
unstable manifold of a saddle point can intersect itself as 
well as the unstable manifold of another saddle point of 
the same saddle cycle. At tr = 0.25 the unstable mani- 
fold of a saddle point intersects the unstable manifold of 
another saddle point of the same period-4 saddle cycle 
(see Fig. 09(b)). Both types of self-intersections appear at 
tr « 0.28 (see Fig. 03(a))- Since our map is non-invert ible, 
a point on the unstable manifold can have two distinct 
preimages so the unstable manifolds are self-intersecting, 
though they are connected. From the definition of unsta- 
ble manifold, we can say that there are infinite number 
of self-intersections (for both types) because the image of 
a self-intersecting point is again a self-intersecting point. 
With increase of the parameter, at tr « 0.28 loops are 



0.02 




-0.04 x -0.01(b) 



FIG. 5: (a) Phase portrait of the attractor at tr = 0.28 with 
TL = -0.3, S L = -0.3, 5r = 1.4 and fi = 0.05. (b) Magnified 
part of the IC marked by the square in (a) displays the loops 
onlC. 



created on the invariant curve (see Fig. 03(a))- Since at 
the points of self-intersections (first type) loops are de- 
veloped so the infinite loop structure is formed on the 
invariant curve (see Fig. 03(b)). The torus formed by this 
invariant curve is a loop torus [20]. As the parameter is 
increased to tr « 0.2855, the eigenvalues of the stable cy- 
cle become complex conjugate. Then the closed invariant 
curve is formed by the union of the unstable manifolds of 
period-4 flip saddle cycle and the points of the period-4 
spiral attractor and flip saddle cycles. Again increas- 
ing tr, the unstable manifolds of the period-4 flip saddle 
cycle nontransversally touch the stable manifolds of it. 
Then the first homoclinic tangency appears at tr « 0.4 
(see Fig. [61(a)). Since the stable eigenvalue of the sad- 
dle cycle is negative, the homoclinic orbit flips around 
the saddle point. Thus the homoclinic tangency occurs 
in both sides of the stable manifold. With further in- 
crease of tr, the stable and unstable manifolds of the 
period-4 flip saddle cycle intersect transversally to form 
a homoclinic tangle (see Fig. EJb)). The existence of ho- 
moclinic tangle gives rise to Smale horseshoe dynamics, 
consequently an infinite number of high periodic orbits 
are created. The torus no longer exists, but the period-4 
spiral attractor and the flip saddle cycle persist. Further 
increasing tr, the second homoclinic tangency occurs at 
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FIG. 6: Phase portraits of the attractors at (a) tr = 0.4 and 
(b) tr = 0.4295 with other parameters same as in Fig. [3] 
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FIG. 7: Attractors are depicted at (a) tr — 0.5 and (b) 
0.53 with other parameters same as in Fig. [3] 



« 0.5 (see Fig. 0(a)). The period-4 spiral attractor 
coexists with a chaotic attractor between the parame- 
ter values rjj; and t£. At « 0.57414 the period-4 
spiral attractor and flip saddle cycles merge and disap- 
pear through border collision fold (saddle-focus) bifurca- 
tion. This scenario of transition to chaos is quite different 
from that in invertible maps since the self-intersections 
of the unstable manifolds and creation of loops on the in- 
variant curve are impossible in invertible maps [l6|, [l9[ . 
There is only one attractor before the first homoclinic 
tangency, namely, the closed invariant curve consisting 
of the periodic orbits. After the appearance of the ho- 
moclinic tangle, the closed invariant curve is destroyed, 
yet the period-4 spiral attractor continues to exist. The 
stable manifold of the saddle point to the left side of 
the border (x = 0) touches the critical curve LC at 
tr « 0.4295 (see Fig. EJb)). Then a contact bifurcation 
[23( | occurs. It causes changes in the basin of attraction 
of the period-4 spiral attractor, simply connected to non- 
connected and unbounded to bounded, since the basin 
boundary is formed by the stable manifolds of the sad- 
dle cycle. A different basin of attraction for the period- 
4 attractor is created. The stable manifolds of each of 
the four saddle points are disconnected from each other. 
The stable manifold of a non-invertible map can inter- 
sect itself when the map is structurally unstable. Ow- 
ing to the presence of the homoclinic tangle, our map 



has already become structurally unstable. So the stable 
manifolds of the period-4 flip saddle cycle form noncon- 
nected (with finite number of simply connected compo- 
nents) and bounded basin of attraction of the period-4 
spiral attractor. To our knowledge, this phenomenon in 
PWS dynamics has not been reported earlier. The stable 
manifolds fold at every intersection with the y-axis Q. 
Other fold points on the stable manifolds are preimages 
of the fold points on the 7/-axis in Z^. As we increase t#, 
the size of the basin of attraction of the period-4 spiral 
attractor diminishes (see Fig. 0(b)). This size tends to 
zero just before the border collision fold bifurcation. 



VI. CONCLUSION 

The study of non-invertible dynamical systems is an 
important and recently flourishing research subject. In 
contrast to the earlier investigations 0, [TBI, Qi| [TtJ , here 
we have chosen the situation where the map (pQ) is non- 
invertible, and the determinant of the Jacobian matrix 
in one side of the bifurcation point is less than unity and 
is greater than unity in the other side. Such a situation 
arises in many physical systems like power electronic cir- 
cuits. In this paper we have illustrated a scenario of tran- 
sition from phase-locked dynamics to chaos in such a non- 
invertible piecewise smooth map. The closed invariant 
curve is formed by the union of the unstable manifolds of 
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a saddle cycle and the points of stable and saddle cycles. 
Since a point can have two distinct preimages, the unsta- 
ble manifolds have structurally stable self-intersections. 
We have found that loops are created on the invariant 
curve although there is no chance of appearance of cusps 
on it. As our map is linear in both sides of the border and 
the critical curve is the x-axis which has no cusp point, 
so cusps do not appear on the invariant curve. Here the 
invariant curve is destroyed through a homoclinic bifur- 
cation. We have shown that contact bifurcation occurs 
when the stable manifold touches the critical curve. This 
bifurcation causes a qualitative change in the structure 
of basin of attraction, from simply connected and un- 
bounded to nonconnected and bounded. At that point, 



the stable manifolds are self-intersecting since the system 
has become structurally unstable after the appearance of 
homoclinic tangle. We have observed the hysteretic tran- 
sitions between mode-locked dynamics to chaotic dynam- 
ics and vice versa and have found the mechanism to be 
quite different from that in invert ible maps. 
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